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A generalization of non-perturbative QED model for high harmonic generation is developed for 
the multi-mode optical field case. By introducing classical-field-dressed quantized Volkov states 
analytically, a formula to calculate HHG for hydrogen-like atom in ultrashort intense laser pulse 
is obtained, which has a simple intuitive interpretation. The dressed state QED model indicates a 
new perspective to understand HHG, for example, the presence of the weak even-order harmonic 
photons, which has been verified by both theoretical analysis and numerical computation. Long 
wavelength approximation and strong field approximation are involved in the development of the 
formalism. 



High order harmonic generation (HHG), that atoms 
subjected to an intense laser field emit more than a 
hundredth-order harmonic photons, is the highest non- 
linear optical phenomenon ever observedpj. It has trig- 
gered a series of experimental and theoretical researches. 
HHG not only offers an approach to the table-top source 
for coherent XUV and soft x-ray generation Q, but also 
to attosecond light pulses 3 . The latter has been used 
to probe ultrafast atomic and molecular dynamics, eg. 
atom ionization and correlation, molecule rearrangement, 
molecule dissociation and so on^. 

HHG has been widely understood as a three-step pro- 
cess — the ionization, propagation and recombination. 
Calculations have been carried out by direct numerical 
integration of the time-dependent Schrodinger equation, 
using base expansion or grid methods [B|. In addition, 
various theoretical models concerning essential physical 
processes in intense laser field and less computer intensive 
are widely studied. For example, the KFR model for the 
scattering^, the ADK model for the ionization^, the 
Lewenstein model for HHG [9], and so on. These mod- 
els are based on the classical treatment of the laser field. 
In these semiclassical models, the HHG spectrum is ob- 
tained by the Fourier transform of the time-dependent 
electric dipole, or the dipole acceleration. 

On the other hand, the observed discrete HHG spec- 
trum indicates that it is an optical-conversion process in 
which q pump photons are combined to produce a single, 
harmonic photon This picture favors the QED ap- 
proach which treats the photon field as a quantized field. 
It is capable of presenting a clearer picture of how the 
photons are being absorbed and emitted, and how the 
energy and momentum transfer between the optical field 
and the atom. A non-perturbative QED scattering the- 
ory of single- mode case has been developed for HHG [13]. 
It provides a new perspective to understand HHG, for 
example, rather than by extracting HHG spectrum from 



the time-dependent atomic wavepacket, it self-contains 
harmonic photon emitting in the theory. Similar phe- 
nomena at much higher intensities can be calculated on 
its framework. It has successfully reproduced the salient 
characteristics of experimental spectra and revealed the 
connection between high harmonic generation and above 
threshold ionization 1( 
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The essential part of QED model for HHG in a intense 
laser field is the transition matrix element [lfj|. 11, T3] 



(i) 
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where ^ is the initial state of the atom-photon system, 
if? f is the final state with a harmonic photon, ^ is a 
quantized-ficld Volkov state, and V (V ) is the interac- 
tion potential. The symbol 5^ represents energy conser- 
vation condition in the non-relativistic case. Eq. ([1]) gives 
a physical picture that the atom in an intense optical field 
can simultaneously absorb more than one photon to tran- 
sit to the quantized-field Volkov states, and finally emit 
a harmonic photon to transit to the final state. 

The quantized-field Volkov states are the pivotal part 
of the QED approach. Classical Volkov statesQ, which 
are the solutions of time-dependent Schrodinger equa- 
tion for a free electron in a classical electromagnetic 
field, are incorporated in the KFR theory to form a 
non-perturbative scattering framework for multiphoton 
ionization (MPI). Later, Guo, Aberg and Drake ob- 
tained analytical solutions for an free electron in an in- 
tense quantized optical field (lR fill ] — the quantized-field 
Volkov states, based on which the QED calculation for 
HHG is developed. 

Most of the calculations for HHG in QED have been 
carried out by assuming the fundamental laser field to 
be single- mode 1(3, U, 12]. The formal scattering the- 
ory, which takes the interaction time as infinity, is always 
used. This is a good approximation to many experiments 
in which the laser pulse's duration is dozens of optical cy- 
cles long so that the laser field can be considered as single- 
mode. However, it still limits the usage of the method in 
certain ways. For example, the effect of the laser pulse's 
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FIG. 1: HHG spectrum of H atom for the laser parame- 
ters: A = 0.8/im, peak intensity / = 3.5 x W 14 W/cm 2 , 
gaussian pulse shape with FWHM = 16/s. The solid 
line is Is — Volkov state — Is transition, the dash line is 
Is — Volkov state — 2p transition, and the dot line is Is — 
Volkov state — 2s transition. The inset shows a part of the 
spectrum from harmonic order 33 to 55. 



shape cannot be studied, the calculated HHG spectrum 
only contains the photon frequencies which are multiples 
of the fundamental frequency, and ultrafast atom-laser 
interactions are hard to describe. 

In the present study, we improve the QED model to the 
multi-mode case by taking the Born-Oppemheimer ap- 
proximation, that the photon mode is a far more rapidly 
varying degree of freedom than the electron, so the stimu- 
lated electron-photon interactions dress the atomic elec- 
tron. Therefore, in the multi-mode case, the interac- 
tion is between the photons and the photon-cloud-dressed 
electron, instead of the naked electron as in the previous 
models. Similar idea has been used in studying laser- 
assisted bremsstrahlung in highly intense laser field 

Two key techniques are used in our theory. Firstly, we 
treat the strongest mode of the laser field as the quan- 
tized mode, which has the same optical cycle as the laser 
pulse. By deriving the classical-field-dressed quantized 
Volkov states, the effects of all the other modes are incor- 
porated as a classical distortion field which dresses the 
quantized Volkov state. Secondly, we make use of the 
Fourier transform of the evolution operator, instead of 
the time-independent formal scattering formalism. This 
allows us to include time in the theory. The inverse op- 
erators are transformed to computable integrals by the 
equations we derive. 

We 
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differ from the previous studies [K 
using coherent states instead of Fock states to 
scribe the laser field. Coherent state |aia 2 ---), 

\a k \ 2 00 n 

where \a k ) — exp ~ J] -r=^\ n ) is the eigenstate of 

n— n ' 

the annihilation operator with dk\ctk} = ak\ajA, 
is a better approximation to the real laser field [171]: 
(a l (t)a 2 (t)---\A\a 1 (t)a 2 (t)---) = £A<(i) = A(t), 
where A(t) — eA(t) is the electromagnetic field vector. 



In this treatment, the absolute value and the phase of 
ctk are in accordance with the amplitude and the phase 
of the kth mode of the laser field. Therefore it provides 
a way to characterize the time-dependence of the laser 

T 

pulse by ctk = y f dre lWkT A{t). The long wavelength 

limit is taken in the present scheme. Walser et al.(l8j 
have considered the multipole contributions to HHG, and 
found such effects be small below an intensity of about 
10 17 W/cm 2 . 

In the non-relativistic case, the total hamiltonian of a 
hydrogen-like atom in the optical field is H = Hq + Vt, 



where Hq 
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is the noninteraction part of the Hamiltonian. U is the 
atomic binding potential; N' and are photon number 
operators of harmonic mode with frequency lo and the 
kth laser mode with frequency w/c, respectively. Vt is the 
electron-photon interaction term Vt = — -^-p-A + j— A . 
Due to the weakness of the harmonic mode, Vt can be 
divided as V T = V + V 
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For simplicity, we denote A as the harmonic mode, A c 
as the quantized laser mode, and A s = ^ A s as all the 

s^c 

other laser modes, so A = A c + A s + A . In general, 
A u = <? ll (e«a ll e lku ' r + c.c), where the polarization vector 
e is given by e = cos(^/2)e x + isin{^/2)l yi where £ is a 
measure of the degree of ellipticity, then £ = and ir/2 
corresponding to linear and circular polarization, respec- 
tively [ld |. The driving optical field in our calculation is 
linearly polarized. 

The non-perturbative approach for HHG is based on 
the fourier transform of the evolution operator U(t,t ) 
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where r = t — t , 0(t) = 1 when r > and Q(t) = 
when t < ; e — > 0+. The operator can be expanded, 
giving the matrix element for the harmonic generation 
process as [13, EH, EH 

^ = (^(7)1^,0)1^(0)) 

= / dze-" T e^ E t+i w ) T x ($/, cv/(T), l'| 
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z — Hq + ie z — H + ie z — Hq + ie 
X \$i,Oi(p), 0') , (5) 
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where the system's initial and final states are decoupled 
states of the electron and the laser field, and the laser 
pulse is confined to the interval [0,T]. <£V/ is the wave 
function of the atomic electron with binding energy Ef ^ ; 

\n ), with n = 0, 1, is the harmonic photon Fock state, 

and \a>ij(t)) = 



is the coherent laser state . 

If <I> is the eigenstate of the atomic hamiltonian and Hq 
is the noninteraction hamiltonian, we get the relations 
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Then the transition matrix TV. becomes 
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Therefore, the two inverse operators have been trans- 
formed into computable integrals. 

To apply this strategy further to tackle Eq. ©, we 
need the knowledge of the eigenstates of H . Due to the 
weakness of V and the strong field approximation(SFA), 
the binding potential U can be neglected for continuum 
electron states hence only the leading term of the in- 



verse operator 



z-H+ie 



has to be kept, thus enabling U 



and V in the denominator in Eq. © to be dropped [13]. 

In the one-mode case, the eigenstates of such a Hamil- 
tonian are the quantized-field Volkov states [lij]: \& pn — 



i In 

Ve 1 YL e^ p+ ^ u "-^- r x Jj(C,V,<t>ty x e~ i ^\n + j) 

j = -n 

r> 2 

with corresponding energy eigenvalues E pn = + 

(n + 1/2 + u p )w. Up = *L A is the ponderomotive en- 
ergy in units of the photon energy of the laser, where 
A, as the limit of g-\fn{g — > 0, ->/n — > oo), is the 



half amplitude of the classical field; £ 
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r\ = ^UpCos^ and <f>£ — arctan(£ Ji tan |). The gen- 
eralized Bessel function Jj{C,,r},(j}^) is 4>z) = 

oo 
m— — oo 

In the multi-mode case, with one specified quan- 
tized mode (A c ), we obtain the classical- field-dressed 
quantized-field Volkov states 

*^(*) = e <eEA ' (t) - r |*^(t))|a'(t)), 

OO 

= V~ 1/2 e ie ^ Asi - t) ' r \a s {t)) ^2 e i|p+ K-jW" 

j=-n c 

x J^CcVcAlYe-^W + j)e- iE ^ . (9) 

As before, the script s represents all the other modes in 
the laser field except the mode distinguished by c. In 
the large-photon- number limit, where a + |a) = a*] a) for 
\a\ — > oo, ty pn has the property 
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Where E s = J2 w A a s\ 2 - Eq. ([TO]) verifies that the 
classical-field-dressed quantized-field Volkov state is the 
approximate eigenstate of the electron-laser subsystem in 
the multi-mode field. Therefore, a relation similar to Eq. 
([6]) can be obtained 
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Combining Eq. J5J) and (fTTj) . and after the integration 
of the energy z, we get 
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where f ( n ) = ^A 2 (n) - £A(n) • (p + e£ A s (n) + - j)k c ) and V'fa) = j^g'e'* ■ A(r 2 ) - %g'e'* ■ (p 
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FIG. 2: One Part of HHG spectrum of H atom for the same 
laser parameters as in Fig. {T}. It only shows HHG from 
Is — Volkov state — Is transition. 

e EA s (r 2 ) + (u;-j)k, 

Eq. (fT2|) indicates that the probability amplitude of 
emitting a harmonic photon is the sum of the probabil- 
ity amplitudes of all the individual events in which the 
electron is ionized at some time and recombines with the 
nucleus to emit the photon at a later time. The picture 
is intuitive, but is not trivial for these probabilities add 
up in an interference way. 

HHG of H atom in an intense ultrashort laser pulse 
has been calculated as a numerical example. The laser 
parameters are: A = 0.8/im, peak intensity I = 3.5 x 
10 u W/cm 2 , gaussian pulse shape with FWHM = 16/s, 
that is approximately 6 optical cycles. The atom is as- 
sumed to be initially in the ground state. Different lines 
in Fig dl]) represent different final state cases, here the Is, 
2p, 2s states respectively. It shows that, among the three 
transitions considered here, Is — Volkov states — Is tran- 
sition dominates the harmonic photon generation. The 
contributions from other transitions are about 10 -1 ~~ 2 
smaller in the plateau. It demonstrates that the SFA, 
where only the recombination to the ground state is taken 
into account, is a good approximation in this case. This 
QED model can then be used to complement other mod- 
els by indicating what bound states should be included 
in the approximation. 

The fine structures of HHG in Fig ([2]) show the gen- 
eration of even-order harmonic photons. The presence 
of even harmonics can be explained as the contribution 
of 5=7 A , which can not be considered as a perturba- 
tion as in the weak light field case. Transition through 
— — p • A changes the atomic state's parity, while tran- 

sition through t^A~ maintains it. Therefore, the argu- 
te 2m ' G 

ment for the absence of even harmonics, that is the se- 
lection rule for the dipole transition, does not hold here. 
In this laser parameter condition, the intensity of the 
even harmonics is much weaker than the odd ones, with 
a factor 1CP 4 . This is in accordance with the experi- 
ments where only odd harmonics are observed. How- 
ever, in the semi-classical model, or numerically solv- 



ing TDSE, the information of even harmonics is often 
buried by the numerical noise or more importantly, by 
other weaker transitions, like Is — Volkov states — 2s and 
Is — Volkov states — 2p shown here. 

In conclusion, a generalization of non-perturbative 
QED model for HHG has been developed to the multi- 
mode case. In our QED model, HHG produced by ul- 
trashort laser pulse can be calculated for hydrogen-like 
atoms. A numerical example has been calculated for H 
atom in intense ultrashort laser pulse and it reveals a 
clear evidence of even harmonics generation. 
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